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i This talk reviews progress in the (semi-) analytic calculations of the thermodynamics of the quark- 

gluon plasma. I shall explain how weak coupling techniques can allow us, through appropriate 
i resummations, to deal with particular non perturbative aspects of the quark-gluon plasma. Then 

I shall show how the exact renormalization group may provide insight into the physics of this 
| multi-scale system. 
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1. Introduction 

There is presently much discussion going on in order to understand the origin of the strongly 
coupled character of the quark-gluon plasma revealed by the RHIC data on heavy ion collisions 
[1]. These data have indeed led to a paradoxical situation. 

On the one hand, as a natural consequence of asymptotic freedom, one expects hadronic matter 
to turn, at high temperature and density, into a gas of quarks and gluons whose free motion is only 
weakly perturbed by their interactions. Weak coupling calculations (based on resummed QCD per- 
turbation theory) reproduce lattice results for temperatures greater than 2.5 to 3 T c [2], and provide 
support to a quasiparticle picture: according to these calculations, the dominant effect of the in- 
teractions is to turn (massless) quarks and gluons into weakly interacting (massive) quasiparticles. 
Recent lattice calculations, that can probe arbitrarily large temperatures, demonstrate the approach 
to the Stefan-Boltzmann limit in a convincing way, in good agreement with weak coupling cal- 
culations [3]. The calculations of the fluctuations of conserved charges (such as baryon number, 
electric charge, strangeness) provide another evidence that the bulk quark behavior resembles that 
a free gas above the deconfinement transition [4]. 

On the other hand, RHIC data do not provide any evidence for ideal gas behavior [1]. Indeed, 
the strong opacity of matter to hard parton propagation (see e.g. [6] for a recent discussion), the 
strong collective elliptic flow, and the low value of the ratio r\/s (where 7] is the shear viscosity and 
s the entropy density) [7], are observations that picture the quark-gluon plasma produced at RHIC 
as a strongly coupled system. 

One may argue that, at RHIC, the quark-gluon plasma spends most of its existence in a region, 
say between T c and ~ 3T C , where the physics is hard and poorly understood. It seems indeed that, 
in this region, the quasiparticle picture breaks down, and genuine non-perturbative effects appear 
in bulk thermodynamics. For instance £ — 3P, the trace of the energy momentum tensor is non 
vanishing, reflecting the QCD scale anomaly. Still, even in that region, explicit calculations reveal 
that the coupling constant is not huge (near T c , g ~ 2, see e.g. [8]). 

At this point, it is worth recalling that non perturbative features may arise in a system from 
the cooperation of many degrees of freedom, or strong classical fields, making the system strongly 
interacting while the elementary coupling strength remains small. In the quark-gluon plasma coex- 
ist degrees of freedom with different wavelengths. Whether these degrees of freedom are weakly 
or strongly coupled depends crucially on their wavelength: short wavelength degrees of freedom 
may be weakly coupled if the coupling constant is small, while very long wavelength modes will 
remain strongly coupled, however small the coupling constant may be. 

These remarks will help us to understand that, while strict perturbation theory breaks down 
even at moderate values of the coupling constant, weak coupling techniques with appropriate re- 
summations can lead to accurate calculations of thermodynamical properties of the quark-gluon 
plasma. This is because short wavelength degrees of freedom dominate the thermodynamics: the 
long wavelength modes carry individually little energy and they have limited phase space, so their 
total contribution to the pressure is small. Short wavelength degrees of freedom can be handled 
with a combination of effective field theory and perturbation theory, as we shall see in the next 
section. 



2 



RG at finite T 



JeanPaul Blaizot 



2. Weak coupling techniques 

2.1 Breakdown of strict perturbation theory 

Much effort has been put into calculating the successive orders of the perturbative expansion 
for the pressure and the series is known now up to order g 6 lng (see [9] and references therein). 
These calculations have revealed that perturbation theory makes sense only for very small values 
of the coupling constant, corresponding to extremely large values of T. For not too small values 
of the coupling, the successive terms in the expansion oscillate wildly and the dependence of the 
results on the renormalization scale keeps increasing order after order (see e.g. [10]), making strict 
perturbation theory inapplicable to estimate the corrections to the ideal quark gluon plasma. Let us 
emphasize that the problem is not specific to QCD, but also occurs in simpler scalar field theories 
(see [10], and also Fig. 3 below). 

This situation is to be contrasted with what happens at zero temperature, where perturbative 
calculations achieve a reasonable accuracy already at the GeV scale. The point is that, as we 
already mentioned, the validity of the weak coupling expansion depends not only on the strength of 
the coupling, but also on the number of active degrees of freedom. At zero temperature, one deals 
most of the time with a very limited number of degrees of freedom (the colliding particles and the 
reaction products), while at finite temperatures, as we shall see shortly, the thermal fluctuations 
alter the infrared behavior in a profound way. 

2.2 The role of thermal fluctuations 

In the quark-gluon plasma, the effect of the interactions at a given scale depends on the mag- 
nitude of the relevant thermal fluctuations at that scale (and in some cases at a different scale as 
well). Let us examine this in detail. In order to avoid being distracted by issues related to gauge 
invariance we shall restrict ourselves to the case of a simple massless scalar field theory, with a 
g 2 <p 4 interaction. The thermal fluctuations are given by the following integral 



When we perform a perturbative calculation, we assume that the "kinetic energy" ~ ((d(f)) 2 ) is 
large compared to the "potential energy" ~ g 2 (0 4 ). Obviously, this comparison depends on the 
strength of the coupling, but also on the typical wavelength, or momentum, of the fluctuations. 
To make things more precise, let us observe that the integral (2.1) is dominated by the largest 
values of k (in the absence of the statistical factor it would be quadratically divergent). One may 
then calculate the integral with an upper cut-off K and refer to the corresponding value as to "the 
contribution of the fluctuations at scale /c", and denote it by {§ 2 )k- In the same spirit, we shall 
approximate the kinetic energy as ((d<p) 2 ) K « K 2 ((f) 2 ) K . Taking furthermore (0 4 )k~(0 2 )k> one 
gets as expansion parameter 



Let us then examine this parameter for several characteristics momenta. The fluctuations that 
dominate the energy density at weak coupling correspond to the plasma particles and have momenta 




(2.1) 
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k~T. For these "hard" fluctuations, 

K~T, (0 2 ) r ~r 2 , y T ~g 2 - (2.3) 

Thus, at this scale, perturbation theory works as well as at zero temperature (with expansion pa- 
rameter ~ g 2 y or rather a = g 2 /4k). 

The next "natural" scale, commonly referred to as the "soft scale", corresponds to K ~ gT. We 

have 

K~gT, ((l> 2 ) g T~gT 2 , 7 g T~g. (2.4) 

In calculating ((f) 2 ) K for K <C T, we have used the approximation « T/k, so that ((f) 2 ) K <^T ~ kT. 
We note that y g r ~ g, so that perturbation theory can still be used to describe the self-interactions 
of the soft modes. However the perturbation theory is now an expansion in powers of g rather than 
g 2 : it is therefore less rapidly convergent. The emergence of this new expansion parameter is the 
origin of odd powers of g in the perturbative expansion of the pressure (such as the plasmon term 

Another phenomenon occurs at the scale gT. While the expansion parameter y gT that controls 
the self-interactions of the soft fluctuations is small, the coupling between the soft modes and 
thermal fluctuations at scale T is not: indeed g 2 (<p 2 )r ~ (gT) 2 . Thus the dynamics of soft modes 
is non-perturbatively renormalized by their coupling to hard modes. This particular coupling is 
encompassed by the so-called hard thermal loops [11]. 

Finally, there is yet another scale, the "ultra-soft scale" K ~ g 2 T, at which perturbation theory 
completely breaks down. At this scale, we have indeed 

K~g 2 T, (<j> 2 ) gT ~g 2 T 2 , Y g 2 T ~\. (2.5) 

Thus the ultra-soft fluctuations remain strongly coupled for arbitrarily small couplings. Of course, 
this situation does not really occur for a scalar field since a mass is generated at scale gT which 
renders the contribution of the fluctuations at the scale g 2 T negligible. However this situation is 
met in QCD for the long wavelength, unscreened, magnetic fluctuations. 

These considerations suggest that the main difficulty with thermal perturbation theory is not so 
much related to the fact that the coupling is not small enough (for the relevant temperatures it is not 
huge, as we have already pointed out), but rather to the interplay of degrees of freedom with various 
wavelengths, possibly involving collective modes. Our main concern here is bulk thermodynamics, 
dominated by hard degrees of freedom. As we shall see in the next subsections, this can be handled 
adequately by weak coupling techniques involving appropriate reorganizations and resummations 
of the perturbative expansion. We shall consider two examples of such approaches. 

2.3 Effective theory- dimensional reduction 

A powerful technique to handle situations where modes at different scales couple is that of 
effective theories. In the present context, it is natural to isolate the mode carrying zero Matsubara 
frequencies, which leads to the so-called dimensional reduction. In its more elaborate version [12], 
this consists in writing a sequence of effective theories, obtained by integrating successively the 
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fluctuations at scale 2tcT and gT. The coefficients of the resulting effective lagrangians can be 
determined perturbatively as a function of the gauge coupling g. 
The QCD partition function can be written as (V is the volume) 

Pqcd = T ~ lnZ QCD , Z QCD = J DAIDA^-S^Q™, (2.6) 

where, to simplify the writing, we ignore the gauge fixing terms. By integrating out the modes with 
momenta ~ 2nT, one can re-write this as 

P Q CD = P E + ^nZ E (2.7) 

where the lagrangian Jt? E in Z E is is that of the effective theory at scale gT, and is of the form 

J% = ^TrF% + Tr[D*, A ] 2 + m 2 TrA 2 + X E (TrA 2 ) 2 + • • • (2.8) 

where the dots denote higher dimension operators. Here the fields depend only on the spatial 
coordinates, Dk = dk — ig E A^ and the various parameters are determined using perturbation theory. 
In leading order, m| ~ g 2 T 2 , g| ~ g 2 T, X E ~ g 4 T, while P E ~ T 4 . 

Integrating the modes at scale gT leaves an effective theory for the modes at scale g 2 T\ 

T - lnZ E =p M + T - lnZ M , Z M = J DA a k e^ 3 *^, Jf M = l -TvF 2 + • • • (2.9) 

In leading order, pm ~ m\T , g 2 ^ ~ g\. This effective theory is equivalent in its leading order to a 3 
dimensional Yang-Mills theory. It gives a non perturbative contribution to the pressure of order g 6 . 

Calculations based on this scheme have been pushed to high order [9], but the determination 
of the order g 6 contribution to the pressure depends on an as yet undetermined 4-loop matching 
coefficient. By adding a parameter to account for this uncalculated contribution, one can match the 
four-dimensional lattice results in a range of temperatures between 2.5T C and 3T C . It is interesting 
to note that better results are obtained (in particular a weaker dependence on the renormalization 
scale) if one refrains from expanding the parameters of the effective theory in terms of the original 
gauge coupling g [13]. 

2.4 Skeleton expansion- 2PI approximations 

There exists other ways to reorganize the perturbative expansion (for a review see [10]). I 
shall briefly recall here an approach based on a ^-derivable two-loop approximation (also called 
2PI formalism) applied to the calculation of the entropy density [14, 2]: 

^ = - r /iw^[ Imlo « D "'- ImnReD ] 

f d 4 k df(co) r , , 

-2tr J (27r)4 J d K T > [imlogS" 1 - ImEReS] , (2.10) 

with the relevant Feynman diagrams displayed in Fig. 1. In contrast to dimensional reduction, 
based on the imaginary time formalism, this approach exploits the real-time information that one 
has on the elementary excitations of the quark-gluon plasma: the spectral information enters the 
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Figure 1: Skeleton (2PI) diagrams that are taken into account at 2-loop order in QCD for the calculation of 
the entropy according to Eq. (2.10). Wiggly, plain, and dotted lines refer respectively to gluons, quarks, and 
ghosts. 



calculation of the entropy through the propagators of gluons (D) and quarks (5), and the respective 
self-energies (Ft and E). It was shown in Refs. [14, 2] that the lattice results for the entropy of the 
gluonic plasma were quite well reproduced for T >3T C (see Fig. 2). 

The formalism used in this calculation of the entropy has been tested [16] in the limit of a large 
number of quark flavors, which can be solved exactly. One then found that the 2PI approximation 
scheme allows for a smooth extrapolation that is accurate up to quite large values of the coupling 
constant. 

This approach has also the virtue of providing a clear physical picture: at high T, the dominant 
effect of the interactions is to turn the original (massless) degrees of freedom, quarks and gluons, 
into massive quasiparticles, with weak residual interactions. As shown by Fig. 2, this picture 
is consistent with lattice calculations for temperatures above 3T C . Note that in the region below 
2.5T C both this approximation and that based on dimensional reduction systematically overshoot 
the lattice data, suggesting that different physics is involved there. 

2.5 Strong coupling and Supersymmetric Yang Mills theories 

A radically different approach consists in assuming that all the degrees of freedom of the 
quark-gluon plasma are strongly coupled, and take this strongly coupled quark-gluon plasma as 
a starting point for further approximations (such as a strong coupling expansion). This is made 
possible by a theoretical breakthrough that allows calculations in (some) strongly coupled gauge 
theories, based on the so-called AdS/CFT duality (for a recent review, see e.g. [17]). One prediction 
of such calculations concerns the entropy S of Jf = 4 supersymmetric Yang-Mills (SYM) theory 
which behaves, in strong coupling, as [18] 

H + I« (3) i^ (2U) 

where X = 2g 2 N c and So is the entropy of the non interacting system. Thus, in the limit of strong 
coupling, X — > oo, the entropy is bounded from below by the value S/So = 3/4, a value that is not 
too distinct from that obtained from QCD lattice calculations for temperatures above 3T C . In fact, 
in this temperature range, the entropy density is about half-way between its weak coupling value 
and its strong coupling value, and there seems to be no compelling reason to favor an interpretation 
of lattice data based on strong coupling (see [19] for a recent discussion). 
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Figure 2: Comparison of the lattice data (grey band) for the entropy of pure gauge SU(3) theory [15] with 
two successive approximations (solide and dashed lines) in the two-loop 2PI formalism (from [14]). 



3. Insights from the exact renormalization group 

As the preceding discussion has shown, the quark-gluon plasma is a multi-scale system. When 
the coupling is small, a clean hierarchy of well separated momentum scales emerges (T,gT,g 2 T). 
This allows us to treat the very high temperature plasma using a combination of effective theories 
and perturbation theory. However, as the temperature decreases, the coupling grows and the various 
scales start to mix, making the accuracy of such calculations worse and worse. In this last part of 
the talk, I shall indicate how the non perturbative (sometimes called exact) renormalization group 
(for a review, see [20]) can be used to handle this situation. The present discussion is limited to the 
case of the scalar field theory: as we have seen, the difficulties of perturbative calculations in QCD 
have much in common with the corresponding calculations for the scalar field; furthermore we have 
for this scalar case results corresponding to rather elaborate solutions of the flow equations. (For a 
recent discussion of the application of the exact renormalization group to QCD see e.g. [21].) 

3.1 The exact renormalization group 

There is some analogy between the effective field theory approach, such as that mentioned 
earlier, and the non perturbative renormalization group: in effective field theory one integrates 
out degrees of freedom above some cut-off; in the renormalization group this integration is done 
smoothly. In a way, the renormalization group builds up a continuous tower of effective theories 
that lie infinitesimally close to each other, and are related by a renormalization group flow equation. 
This picture is independent of the value of the coupling, so that the renormalization group provides 
a smooth extrapolation from the regime of weak coupling, characterized by a clean separation of 
scales, towards the strong coupling regime where all scales get mixed. In a recent study [22], it 
has been shown in the case of a scalar 4 theory with 0(N) symmetry, that such a technique can 
provide a smooth extrapolation to strong coupling, which turns out to be similar to that of a simple 
(2-loop) 2PI approximation. 
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Figure 3: The pressure as a function of the coupling constant 

In practice, the exact renormalization group builds a family of theories indexed by a momen- 
tum scale parameter fc, such that fluctuations are smoothly taken into account as K is lowered from 
the microscopic scale A down to 0. This is achieved by adding to the original Euclidean action 
S a term of the form AS K [q>] = \$ q Rk(<1)<P(<i)<P(—<i), where the cut-off function R K (q) is chosen 
so that: i) R K (q) ~ K 2 for q < K, which effectively suppresses the modes (p(q < k), and ii) R K (q) 
vanishes for q>K, leaving the modes (p(q> k) unaffected. One can write for T K [<p] an exact flow 
equation [23]: 

d K TM = \J ^d K R K (q)G K (q-q;<t>), (3.1) 

where G K [<j>] is the full propagator in the presence of the background field 0: 

G~ 1 [(l>]=r^[<l>}+R K , (3.2) 

with r^[0] the second functional derivative of T K [(j>] w.r.t. (j). The initial conditions of the flow 
equation (3.1) are specified at the microscopic scale K = A where fluctuations are frozen by AS K , 
so that r^A^] ~ S[(j)]. The effective action of the original scalar field theory is obtained as the 
solution of (3.1) for K — > where R K (q 2 ) vanishes. 

Differentiating Eq. (3.1) m times with respect to yields the flow equation for the m-point 
function [q\ , . . . , q m \ (j)]. Thus for instance, the flow equation for r^ 2 ) reads: 

^ifW) = J ^L d d K R K {q)G 2 K {qA) 

x ) (P, q, ~P - q\<l>)G K (q + p, </>)ri 3) (-p, p + q, -q; «/>) - (p, -p, q, -q; <j>) j . 

(3.3) 
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Note that the flow equation for \q\ , . . . , q m \ 0] involves r[ m+1 ^ and r[ m+2 \ leading to an infinite 
hierarchy. Approximations need to be done to close this hierarchy. We shall present results obtained 
with two distinct approximations. 

3.2 Approximations 

The simplest approximation is the so-called local potential approximation (LPA). It leads to a 
closed flow equation for the effective potential V K . This follows from that of the effective action 
F K , Eq. (3.1), when restricted to constant 0. It reads 

Kd K V K (<j>) = l -j ^Kd K R K (q)G K (q,(j>), (3.4) 

where 

G-\q,<t>)=r ( ?(q,$)+R K ( q ). (3.5) 

In the local potential approximation one assumes that (q, 0) = q 2 + d 2 V/d(j) 2 , so that, with this 
ansatz, the equation for V is indeed a closed equation. 

A more refined approximation [24] (BMW), consists in neglecting the ^-dependence of the 
m-point function in the right hand side of the flow equations, while keeping the full dependence on 
the external momenta pi. For instance, in Eq. (3.3), the approximation amounts to the replacements 
T^\-p,p + q, -<?;0) — > r^ 3) (-p, 0; and ^ 4) (p, -p, ^, (/)) — ► ^k\p, -p,O,O;0). The 
hierarchy is then closed by observing that r[ m+1 ^ {p\ , . . . , p mi 0, 0) = d^F^ {p\ , . . . , p mi 0). This 
approximation has been recently tested on the critical 0{N) model and it was found to yield excel- 
lent values for the critical exponents [25]. 

3.3 Results 

We now present briefly results obtained for the pressure of a massless scalar field theory with 
a g 2 4 coupling [26]. 

Fig. 3 displays the pressure as the function of the coupling constant at the scale g(2nT). 
The various diverging curves indicate the results of perturbative calculations, up to order g 6 (for 
recent high order calculations of the thermodynamics of the scalar field see [27, 28]). The blue 
line indicate the result of a 2PI calculation based on a 2-loop skeleton. The black curve indicate 
the result of the LPA with a Litim regulator, while the blue dots represents the results of the full 
BMW calculation with an exponential regulator. In contrast to the perturbative calculation, the 
two calculations based on the renormalization group show a remarkable stability, and a smooth 
extrapolation towards strong coupling. 

This stability is even more obvious on the next plot, Fig. 4. Now the pressure is plotted as 
a function of the thermal mass. This has the advantage of eliminating all scheme dependence 
[22] (scheme dependence is present when the pressure is plotted as a function of the coupling 
constant). Note the excellent agreement between the 2PI calculation and the BMW one. Note also 
the agreement between the LPA and BMW up to large values of the coupling constant. This plot 
was prepared with the thermal mass calculated up to order g 4 , and accordingly the perturbative 
results for the pressure are reported here only up to this order. 
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Figure 4: The pressure as a function of the thermal mass 



The stability of these results as one moves towards larger values of the coupling constant, in 
marked contrast to those of strict perturbation theory, suggests that further corrections to the present 
calculations (within the renormalization group or the 2PI formalism) may indeed be very small. 

Acknowledgements. I thank Andreas Ipp for his help in preparing Figs. 3 and 4. I am also grateful 
to the Graduate School of Art and Science of the University of Tokyo for support and hospitality at 
the time of the write-up of this contribution. 
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